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PROCEDE DE COMMUNICATIONS RADIOMOBILES AMRT ITERATIF 

DESCRIPTION 

5 Domaine i:echn±que 

La presente invention a pour objet un procede de 
coiTununications radiomobiles AMRT ("acces multiple a 
repartition en temps", en anglais TDMA pour "Time 
Division Multiple Access") iteratif. 
10 Elle peut s'appliquer au systeme radiomobile 

europeen GSM [1] et son evolution GSM++ ou au systeme 
americain D-AMPS. 

E^a^ de la t:echn±qpie ant:er±eujre 

15 Les systemes AMRT permettent de repartir dans le 

temps plusieurs utilisateurs communiquant a travers un 
meme canal radioelectrique [1, 3] • On separe chaque 
utilisateur des autres en lui assignant periodiquement 
une tranche temporelle (ou intervalle, ou creneau 

20 temporel, ou "timeslot" en anglais) pour lui permettre 
d'emettre et une autre tranche temporelle pour lui 
permettre de recevoir. Pour se premunir des 
imperfections en terme de synchronisation temporelle 
entre emetteurs et de 1' interference entre tranches 

25 temporelles consecutives produite par la pluralite de 
trajets, des intervalles de garde sont prevus de part 
et d' autre de chaque tranche temporelle. 

La figure 1 annexee montre ainsi des tranches 
temporelles TAi, TA2, TAn allouees a N 

30 utilisateurs, avec des intervalles de garde IG. 
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Le signal radioelectrique emis durant une tranche 
temporelle est obtenu par transposition en frequence du 
signal en bande de base equivalent. Le signal en bande 
de base est le resultat du filtrage du train de donnees 
a emettre par un filtre dit d' emission. Habituellement , 
le train de donnees emis est compose de deux sequences 
de symboles de donnees separees en temps par une 
sequence de symboles connue du recepteur, dite sequence 
de reference. Les deux sequences de symboles de donnees 
peuvent eventuellement provenir du codage et de 
1' entrelacement des informations a emettre. 

Afin de permettre une estimation precise du canal 
au niveau du recepteur, la sequence de reference 
possede souvent de bonnes proprietes de correlation 
temporelle. A titre d'exemple, les sequences de 
reference utilisees dans le cadre du systeme GSM^ 
appelees sequences CAZAC (contraction de "Constant 
Amplitude - Zero Autocorrelation" en anglais) [1, 4, 
5] f sont des sequences a composantes prises dans 
1' alphabet bipolaire {-1^ 1} et possedant une fonction 
d' autocorrelation circulaire nulle partout sauf a 
1' origine. 

Le canal radiomobile emprunte lors d'une 
communication entre un emetteur et un recepteur est 
generalement de type multi-tra jets avec evanouissements 
rapides dits de IIAYLEIGH. L' existence de plusieurs 
trajets vient de ce que 1' onde radioelectrique se 
propage selon plusieurs chemins entre le lieu 
d' emission et le lieu de reception. 

Le signal re?u est alors la somme de plusieurs 
repliques plus ou moins retardees et plus ou moins 
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alterees en phase et en amplitude. Pour obtenir une 
restitution fiable de 1' information transportee par ces 
repliques, le recepteur opere un filtrage adapte au 
filtre d' emission et au canal et il combine de maniere 
optimale les contributions energetiques de toutes les 
repliques du signal emis. Le signal en sortie du filtre 
adapte est echantillonne au rythme des symboles et il 
est blanchi au moyen d'un filtre discret dit 
blahchisseur • 

Les echantillons en sortie du filtre blanchisseur 
fournissent une version filtree et bruitee du train de 
donnees emis. Le filtre discret associe a cette version 
filtree, egalement appele canal discret, presente une 
reponse impulsionnelle finie qui varie d'un echantillon 
a 1' autre. II caracterise d'une maniere indirecte la 
realisation du canal multi-tra jets durant la tranche 
temporelle de reception cbrrespondante . 

Les echantillons en sortie du filtre blanchisseur 
correspondant a la sequence de reference sont utilises 
pour estimer le canal discret [1, 4, 5]. Cette 
estimation du canal discret permet d'egaliser. les 
echantillons restants permettant eventuellement de 
detecter les deux sequences de donnees emises. 
L'egaliseur habituellement utilise est connu sous le 
nom d'egaliseur de VITERBI avec decisions souples [2]. 
Comme son nom I'indique, cet egaliseur fournit, grace k 
une application de I'algorithme de VITERBI avec sorties 
ponderees ("Soft-Output Viterbi Algorithm" ou SOVA en 
anglais) [2], des decisions souples sur tous les 
symboles de donnees emis. Ces sorties souples sont 
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eventuellement desentrelacees et decodees pour detecter 
les informations envoyees . 

La figure 2 annexee illustre ces operations 
d ' emission/reception. L'emetteur E comprend une source 
de donnees 10, un codeur/entrelaceur 12 delivrant des 
symboles ak, un modulateur 14, Le canal multi-tra jets 
est symbolise par le bloc 20. Le recepteur R comprend 
un demodulateur 30 (filtre adapte/f litre blanchisseur) 
delivrant des echantillons Rjc, un estimateur de canal 
discret 32, un egaliseur 34, un desentrelaceur/decodeur 
36 et finalement un destinataire 38. 

Le canal discret vu a la sortie du filtre 
blanchisseur peut varier de maniere significative d'une 
tranche temporelle a la suivante. Cette variation est 
due principalement au changement des conditions de 
propagation entre I'emetteur et le recepteur et a la 
stabilite en frequence au niveau du recepteur. 

Les conditions de propagation ont une influence 
directe sur le canal multi-tra jets observe. Elles 
changent en raison d'une modification de 
1' environnement ou d'un deplacement de I'emetteur et/ou 
du recepteur. Elles entrainent une variation temporelle 
du canal discret, aussi bien entre tranches temporelles 
successives qu'a I'interieur d'une meme tranche. 

La variation du canal discret entre deux tranches 
temporelles successives allouees au meme utilisateur 
est d' autant plus grande que I'intervalle temporel 
entre ces tranches est grand. Cette variation est 
accentuee, meme au niveau d'une seule tranche 
temporelle, en raison d'une augmentation de la 
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frequence porteuse ou de la vitesse de I'emetteur ou/et 
du recepteur . 

Dans la pratique, la variation du canal discret 
entre deux tranches temporelles est suf f isamment grande 
pour empecher toute estimation adaptative de celui-ci, 
L' estimation de canal au niveau d' une tranche 
temporelle doit alors reposer sur les seuls 
echantillons de la sequence de reference 
correspondante . 

L' utilisation de cette sequence de reference 
fournit une estimation non biaisee du canal discret au 
milieu de la tranche temporelle. A I'instar du systeme 
GSM, on fait generalement I'hypothese d'un canal 
invariant au cours d'une tranche temporelle regue. Dans 
ce cas pr6cis^ 1' estimation obtenue au milieu de la 
trame peut servir sans degradation a 1' egalisation du 
reste de la trame. Mais, le deplacement a grande 
Vitesse de certains terminaux et la demande soutenue en 
services a hauts debits operant a des frequences 
radioelectriques de plus en plus grandes, rendent cette 
hypoth^se de moins en moins justifiee. Le canal discret 
peut en effet subir des variations signif icatives entre 
le debut et la fin d'une tranche temporelle donnee. 
L'ecart entre ce canal discret reel et son estimation 
se creuse de plus en plus lorsqu'on s'eloigne de la 
sequence de reference. Ceci peut entralner une 
degradation importante et irreversible de la qualite de 
la reception et/ou des performances du systeme AMRT. 

Une solution [13] a deja ete proposee pour 
resoudre ce probleme de variation rapide du canal dans 
le cas des canaux radioelectriques lies aux 
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radiocoininunications par satellite. Cependant, cette 
solution ne traite que le cas d' un canal avec un seul 
trajet et suppose que les symboles^ de reference ou de 
donnees, ont une energie transmise commune, Elle ne 
peut done pas etre appliquee aux canaux multi-tra jets 
rencontres dans les communications radiomobiles 
terrestres . 

Par ailleurs, a I'instar de 1' extension EDGE du 
systeme GSM, des debits importants peuvent etre 
obtenus, a bande f requentielle inchangee, par la mise 
en oeuvre de modulations utilisant un plus grand nombre 
d'etats. Le systeme devient alors tres sensible a la 
qualite de 1' estimation de canal, et il peut subir des 
degradations importantes meme a de tres faibles 
vitesses de deplacement, si la puissance et/ou la 
longueur de la sequence de reference ne sont pas 
augmentees * 

Le but de la presente invention justement est de 
remedier a ces inconvenients . 

Expose de 1 ' invention 

La presente invention a pour but d'accroitre les 
performances des systemes AMRT en reduisant le rapport 
entre la puissance regue et 1' interference d'acces 
multiple pour une qualite de reception donnee ou, en 
d'autres termes, d'ameliorer la qualite de reception a 
puissance de reception constante. Cette amelioration de 
la qualite permet, entre autres, d'augmenter la 
capacite et la couverture des systemes AMRT, tels que 
le systeme GSM et surtout son extension EDGE pour 
offrir des services a hauts debits. Cette amelioration 
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est obtenue notamment grace a une optimisation de 
1' estimation de canal au niveau du recepteur. Cette 
optimisation permet de contrecarrer toute degradation 
en performances engendree par une variation rapide du 
canal entre le debut et la fin d'une meme tranche 
temporelle. Elle permet egalement de reduire 
signif icativement les degradations engendrees par une 
augmentation du nombre d'etats de la modulation et ceci 
sans recours a une augmentation de la longueur ou de la 
puissance de la sequence de reference. 

Un autre but de 1' invention est de reduire le prix 
de revient des terminaux en les rendant, a qualite de 
reception inchangee^ moins sensibles a 1 ' instabilite en 
frequence d'un oscillateur local bon marche permettant, 
entre autres, la transposition du signal regu en bande 
de base. 

L' invention permet encore de reduire, a qualite de 
reception egale, le nombre et/ou la puissance des 
symboles de reference transmis durant chaque tranche 
temporelle. Ce but est atteint par la prise en compte, 
dans 1' estimation du canal discret, des symboles de 
donnees qui sont en general plus nombreux que les 
symboles de reference. Ce but est egalement atteint par 
la prise en compte de maniere optimale des echantillons 
correspondant aux symboles de reference et/ou de 
donnees d'un nombre arbitraire de tranches temporelles 
allouees ^ 1' utilisateur . Cette prise en compte est 
effective meme en cas de variations importantes du 
canal entre deux tranches successives appartenant au 
meme utilisateur. Ce but est enfin atteint en prenant 
en compte de maniere optimale les echantillons 
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provenant du canal de signalisation commun et 
d'eventuels utilisateurs sur la liaison descendants . 

Selon 1' invention, les symboles de reference sont 
repartis en groupes sur toute la duree de la tranche 
temporelle allouee, ce qui garantit une meilleure 
poursuite de la variation du canal et done une 
robustesse accrue vis a vis des grandes vitesses de 
deplacement et de 1 ' instabilite en frequence de 
1' oscillateur local. Les sequences de reference peuvent 
etre choisies independamment les unes des autres, mais 
elles doivent tou jours presenter de bonnes proprietes 
d' autocorrelation. Elles peuvent etre derivees de 
sequences de type CAZAC de longueur reduite. 

Avec 1' invention, I'energie transmise par symbole 
peut varier d'un symbole a 1' autre sans que cela pose 
des difficultes pour 1' estimation du canal, laquelle 
s'effectuera tou jours de maniere optimale. 

Selon 1' invention, le traitement d'une tranche 
temporelle se fait par bloc a chaque fois que les 
echantillons correspondant a cette tranche et 
eventuellement a ceux d'autres tranches (provenant du 
meme utilisateur, d'autres utilisateurs ou du canal de 
signalisation commun) sont disponibles pour effectuer 
1' estimation de canal. 

A I'instar d'un recepteur AMRT classique (tel que 
celui du systeme GSM) on commence tou jours par une 
estimation grossiere du canal discret en utilisant 
uniquement des echantillons dependant des seals 
symboles de reference de la tranche temporelle a 
traiter. Cette estimation grossiere permet de 
caracteriser, au niveau de chaque symbole de reference 
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ou de donnee du bloc a traiter, symbole par symbole, 
1' evolution de la phase et de 1' amplitude de tous les 
coefficients (egalement appeles "trajets" par abus de 
langage) du canal discret retenus pour la realisation 
5 de 1'' egalisation. 

L'une des caracteristiques de 1' invention est 
1 ' utilisation d'un treillis pour representer le canal 
multitrajets . Cette caracteristique reduit la 
complexite du recepteur. En effet, si I'on prend comme 

10 exemple des sequences de symboles de donnees comprenant 
n symboles, sans codage, avec une modulation de phase a 
deux etats (MDP2) , chaque bit emis peut prendre I'une 
ou 1' autre de deux valeurs, Comme ces valeurs sont 
equiprobables, cela signifie qu'il existe 2" sequences 

15 de symboles de donnees possibles. La complexite d'un 
calcul traditionnel utilisant les probabilites 
partielles pour toutes les sequences possibles crolt 
done exponentiellement (2") avec le nombre de symboles 
par sequence. 

20 L ' utilisation, conformement a 1' invention, d'un 

treillis pour representer le canal reduit 
considerablement la complexite des calculs, car ceux-ci 
vont se subdiviser en sections, chaque section donnant 
lieu a un certain type de calcul, tou jours le meme pour 

2 5 chaque section. La complexite augmente done comme le 
nombre de .sections, lequel croit comme le nombre n. La 
complexite n'est done plus exponentielle mais lineaire. 

Les etiquettes des branches des sections du 
treillis changent au gre de la variation temporelle du 

30 canal discret au niveau de chacune des tranches 
temporelles constitutives du bloc. Par ailleurs, un 
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algorithme different de 1' algorithme de VITERBI avec 
sorties souples (SOVA) [2] est utilise pour une 
caracterisation plus precise des probabilites des 
branches du treillis conditionnellement aux 
echantillons req:us. Dans le cas d'une mise en oeuvre 
optimale de 1' estimateur de canal, un algorithme dit de 
BAHL ou BCJR (abreviation de "Bahl Cocke Jelinek 
Raviv") [6] peut etre utilise pour fournir ces 
probabilites conditionnelles . D'autres algorithmes [10] 
plus simples, provenant generalement d'une 
linearisation de cet algorithme, peuvent etre utilises, 
avec cependant une legere degradation en performance, 
lis permettent de reduire considerablement la 
complexite de 1' estimateur de canal. 

Dans le cas d' un recepteur classique, les seules 
sorties ponderees fournies par 1' algorithme SOVA sont 
directement exploitees apres desentrelacement pour 
decoder les informations transmises. Dans le cas du 
recepteur de 1' invention, les probabilites 
conditionnelles carac]terisant indirectement les 
symboles de donnees sont exploitees, en plus des 
symboles de reference, pour fournir une estimation de 
meilleure qualite du canal discrete 

L' estimation amelioree du canal discret obtenue a 
la fin d'une iteration donnee peut etre utilisee pour 
augmenter la qualite et la fiabilite des probabilites 
conditionnelles des branches du treillis caracterisant 
le canal discret. Ces probabilites conditionnelles 
ameliorees peuvent, a leur tour, etre utilisees, 
conjointement avec les symboles de reference, pour 
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apporter une amelioration supplementaire a 1' estimation 
du canal discret. 

En theorie, une estimation iterative ne conduit a 
une estimation optimale du canal discret qu'apres un 
nombre infini d* iterations . En pratique cependant, 
quelques iterations suffisent sans grande degradation 
des performances par rapport a une estimation optimale. 

Une fois 1' estimation du canal discret accomplie, 
le treillis correspondant peut etre utilise soit par 
I'algorithme optimal de BAHL (ou tout autre 
simplification de celui-ci) , soit par I'algorithme SOVA 
pour fournir des sorties ponderees ameliorees au 
decodeur . 

Le caractere optimal du recepteur selon 
1' invention est lie a la qualite de 1' estimation du 
canal discret vu durant une tranche temporelle donnee. 
Ce caractere optimal repose sur 1' utilisation d' uri 
algorithme iteratif dit SAGE (abreviation de "Space- 
Alternating Generalized Expectation-Maximization 
Algorithm") [7] pour trouver la realisation de canal la 
plus vraisemblable conditionnellement aux echantillons 
req:us du bloc a traiter. L' algorithme SAGE est une 
extension de 1' algorithme EM (contraction de 
"Expectation-Maximization") [8, 9] qui permet de se 
debarrasser du probleme de couplage entre coefficients 
du canal discret lors de leur estimation. 

L* estimation du canal discret repose egalement sur 
la decomposition de chaque trajet du canal discret a 
estimer selon un algorithme d' expansion dit de 
KARHUNEN-LOEVE [3]. Cette decomposition fournit, d'^une 
part, une caracterisation souple des variations 
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temporelles de chacun des trajets discrets et, d'autre 
part, s'integre naturellement dans 1' algorithme SAGE. 
Elle permet egalement de traiter le cas ou I'energie 
emise par symbole varie d'un symbole a 1' autre. 

L' estimation du canal discret repose enfin sur 
1' utilisation de 1' algorithme de BAHL, qui fournit a 
chaque iteration de 1' algorithme SAGE les probabilites 
des branches du treillis caracterisant le canal 
discret, conditionnellement aux echantillons reg:us et a 
1' estimation de canal fournie a 1' iteration precedente. 

De fagron precise, 1' invention a done pour objet un 
precede de communications radiomobiles du type a acces 
multiple a repartition en temps, dans lequel on alloue 
a plusieurs utilisateurs des tranches temporelles pour 
1' emission et la reception de signaux radioelectriques, 
et dans lequel : 

- a 1' emission, on constitue, dans chaque tranche 
temporelle allouee, des symboles de donnees et 
des symboles de reference, on emet un signal 
radioelectrique contenant ces symboles de 
donnees et ces symboles de reference, ce signal 
empruntant un canal radiomobile multi-tra j ets, 

- a la reception, on effectue un filtrage adapte 
du signal regu, on echantillonne le signal 
filtre au rythme des symboles, on effectue une 
estimation du canal radiomobile, on traite les 
echantillons en tenant compte de cette 
estimation du canal et I'on restitue les 
symboles de donnees propres a la tranche 
allouee. 
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ce precede etant caracterise en ce que : 

- a 1' emission, pour chaque tranche temporelle 
allouee, on constitue plusieurs sequences de 
reference possedant des proprietes de correlation 
temporelle appropriees et I'on repartit ces 
sequences dans chaque tranche temporelle, 

— a la reception : 

i) pour estimer le canal radiomobile, on 
represente ce canal par un treillis 
comprenant des branches symbolisant les 
transitions entre deux etats successifs du 
canal et caracterisant les sequences 
possibles a la sortie du canal, 
ii) on calcule les probabilites de ces branches 
en mettant en oeuvre un algorithme iteratif 
dans lequel on calcule, par un nombre fini 
(D+1) d* iterations et selon le critere de 
probabilite maximum a posteriori, une 

estimation d'une representation {g^^^^I^Jo des 
L trajets de ce canal, 1' iteration de rang 
(d+1) (3d+l) permettant d'obtenir une 

reestimation {o^^'^^^^fe^l a partir de 

1' estimation {G^^''^}^rQ obtenue par 1' iteration 
precedente (3d) , et, k chaque iteration 
(3d+l), on calcule en L etapes les vecteurs 
|G^(d+i)j^ ^=0, 1, L-1 a partir des 

vecteurs {g'"^'^*^^ et {G-<^>t, calcules a 
I'etape precedente, 1* estimation (c'^^^j^Jo des 
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L trajets du canal etant finalement obtenue 
a partir de 1 ■ estimation de sa 

representation {g^^^^}^^q . 



5 Breve descrip1:lon des dessins 

- la figure 1, deja decrite, montre la repartition 
des tranches temporelles allouees a plusieurs 
utilisateurs ; 

- la figure 2, deja decrite, est un schema 
10 illustrant 1' emission et la reception dans la 

technique AMRT ; 

- la figure 3 montre les variations du taux 
d'erreur binaire (TEB) en fonction de I'etalement 
Dopp 1 e r no rma 1 i s e ( BdTs ) ; 

15 - la figure 4 montre les variations du taux 

d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour des modulations MDP4 et MDP8, dans une 

premiere configuration et pour un canal discret a 
deux trajets ; 

20 - la figure 5 montre les variations du taux 

d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour des modulations MDP4 et MDP8, dans une 

deuxieme configuration avec un canal discret a 
deux trajets ; 

25 - la figure 6 montre les variations du taux 

d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour une modulation MDP4 utilisant la 

deuxieme configuration et un canal discret a 
quatre trajets ; 
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- la figure 7 montre les variations du taux 
d' erreur binaire (TEB) en fonction du rapport 
Eg/N pour una modulation MDP4 utilisant un canal 
discret avec un, deux et quatre trajets ; 

5 - la figure 8 est une representation du canal 

discret vu a la sortie du filtre blanchisseur, 
sous forme de registre a decalage ; 

- la figure 9 illustre 1' estimation des 
coefficients du filtre du canal discret par 

10 correlation des echantillons regus avec une 

sequence CAZAC et ses decalees circulaires ; 

- la figure 10 montre la repartition des valeurs 
propres les plus importantes de la matrice de 
covariance pour trois etalements Doppler 

15 normalises ; 

- la figure 11 illustre une representation en 
treillis du canal discret normalise vu a la 
sortie du filtre blanchisseur pour un canal a 
deux trajets et une modulation a deux etats et 

20 une longueur de tranche temporelle egale a 

quatre ; 

- la figure 12 montre la chronologie des 
traitements effectues a la reception pour une 
tranche temporelle req:ue ; 

25 - la figure 13 montre la chronologie des 

traitements effectues dans une iteration 
elementaire de I'algorithme d' estimation de 
canal. 
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Descrxp-blon de modes partlcullers de mlse en oeuvre 

Les performances du precede selon 1' invention ont 
ete comparees a celles d'une variants du precede 
classique utilise dans le cadre du systeme GSM [l],Le 
5 recepteur classique du systeme GSM estime d'abord le 
canal discret au milieu de la tranche temporelle a 
traiter^ par correlation des echantillons re<?us 
dependant uniquement de la sequence de reference avec 
la sequence CAZAC correspondante . II egalise ensuite 

10 les autres echantillons rectus au moyen de I'algorithme 
SOVA afin de decoder les symboles d' information emis. 
Afin d'evaluer equitablement 1' amelioration intrinseque 
de performances apportee par 1' invention par rapport au 
recepteur classique du GSM, I'algorithme sous-optimal 

15 SOVA a ete remplace par I'algorithme optimal de BAHL 
utilise dans 1' invention. 

La comparaison est obtenue a travers 1' evolution 
du taux d'erreur binaire (TEB) brut (sans prise en 
compte d'un eventuel codage correcteur ou detecteur 

20 d'erreurs), pour des symboles de reference aussi bien 
groupes au milieu de la tranche temporelle a traiter (a 
la maniere du systeme GSM) que repartis en petits 
groupes sur toute la duree de cette tranche 
(conf ormement a 1 ' invention) , afin de garantir une 
25 meilleure poursuite de la variation temporelle du 
canal. 

La figure 3 represents ainsi 1' evolution du TEB en 
fonction de I'etalement Doppler normalise Bj^T^ , produit 
de I'etalement Doppler et de la duree d'un symbole 

30 Ts , qui caracterise I'intensite de variation du canal 
due au deplacement de I'emetteur et/ou du recepteur. 
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L'etalement Doppler est ^ la fois proportionnel a la 
Vitesse du terminal mobile et k la frequence de la 
porteuse. L' instabilite en frequence de 1' oscillateur 
local, qui n'est pas prise en compte ici, conduit a des 
degradations suppl6mentaires en performances. A 
I'instar du systeme GSM, le comportement du TEB a ete 
evalue pour des tranches temporelles composees de N=14 8 
symboles chacune. La modulation retenue est la 
modulation MDP4 (modulation a deplacement de phase a 4 
etats) . Le canal discret considere est compose de deux 
trajets (L=2) de meme puissance moyenne. Le rapport 
Eg/M entre I'energie moyenne regue par symbole et le 

niveau de bruit est de 10 dB. Trois configurations ont 
ete considerees : 

■ Une sequence de reference centrale de 26 symboles 
bipolaires prise au hasard parmi les 8 sequences 
utilisees dans le systeme GSM. Ces sequences sont 
des extensions cycliques des 8 sequences CAZAC de 
longueur 16 [1, 11] . 

■ Une sequence de reference centrale, {-1, +1, +1, 
+ 1, -1}/ de 5 symboles bipolaires obtenue par 
extension cyclique (a un signe et un decalage 
pres) de 1' unique sequence CAZAC bipolaire de 
longueur 4 a savoir {-1, +1, +1, +1}. 

■ Cinq groupes de sequences de reference, repartis 
sur toute la tranche temporelle a traiter, 
identiques a la sequence de reference centrale de 
5 symboles bipolaires precedente (25 symboles de 
reference au total) . 
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Sur la figure 3, les courbes representees 
correspondent aux cas suivants : 

41 : recepteur classique, premiere configuration, 

42 : recepteur classique, deuxieme configuration, 

43 : recepteur de 1' invention, deuxieme 
configuration, 

44 : recepteur classique, troisidme 
configuration, 

45 : recepteur de 1' invention, troisieme 
configuration. 

Ces resultats montrent que 1' utilisation de 
sequences groupees au centre de chaque tranche 
temporelle (premiere et deuxieme configurations) ne 
permet d'aller au-dela d'un etalement Doppler normalise 
BdTs=1/2500 sans degradation significative des 
performances. A 1' inverse, 1' utilisation de groupes de 
petites sequences reparties (troisieme configuration) 
permet d'atteindre, sans degradation perceptible, des 
etalements Doppler normalises proches de BdTs=1/100. 
Cette amelioration significative de la robust esse du 
recepteur vis-a-vis de I'effet Doppler, et 
eventuellement de 1 ' instabilite en frequence de 
1' oscillateur local, est due a la repartition des 
sequences de reference qui permet de suivre de maniere 
plus precise toute evolution temporelle du canal. 

Le tableau 1 caracterise de maniere chiffree 
1' amelioration en performances apportee par 1' invention 
vis-a-vis de I'algorithme classique ameliore du systeme 
GSM. Dans ce tableau R.C, signifie "recepteur 
classique" et R.I. recepteur selon 1' invention. Le 
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rapport Eg/NO entre I'energie moyenne regue par symbole 
et le niveau de bruit est de 10 dB. Dans le cas d' une 
sequence de reference centrale longue (premiere 
configuration), 1' amelioration apportee par 1' invention 
5 est negligeable. 



Etalement 

DOPPLER 

normalise 


Pi 
con! 


•emigre 
Iguration 


Deuxidme 
configuration 


Troisifeme 
configuration 


R.C. 


R.I. 


R.C. 


R.I. 


R.C. 


R.I. 


1/2500 


0,00 
99 


0,00855 


0,01584 


0,0089 


0.0105 


0,0077 


1/100 


X 


X 


X 


X 


0,0143 


0,0083 



Tableau 1 

10 

La sequence de reference, dans ce cas , est 
suf f isamment longue pour garantir, sans I'aide des 
sequences de donnees inconnues, une bonne qualite 
d' estimation de canal. Cependant, dans le cas d'une 

15 seule sequence de reference courte (deuxieme 
configuration), le recepteur selon 1' invention permet 
d'obtenir non seulement de meilleures performances que 
le recepteur classique, mais aussi des performances 
comparables a ceux, obtenus avec la sequence de 

20 reference longue de 26 symboles de reference. Par 
consequent, le recepteur selon 1' invention est rendu, 
d'une certaine maniere, insensible a une reduction meme 
importante de la taille de la sequence de reference . Ce 
resultat s'explique par la prise en compte, dans 

25 1' estimation de canal, des symboles des sequences de 
donnees qui viennent cooperer avec les symboles de 
reference. Dans le cas de cinq groupes de sequences de 
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reference courtes (troisieme configuration) , le 
recepteur selon 1' invention permet de maintenir le taux 
d'erreur binaire presque inchange en faisant varier 
I'etalement Doppler normalise entre 1/2500 et 1/100. 

5 

Les figures 4, 5 et 6 representent le taux 
d'erreur binaire en fonction du rapport 'E^ / Nq entre 
I'energie moyenne regue par symbole et le niveau du 
bruit pour un etalement Doppler normalise BdTs=1/2500 

10 peu contraignant , 

Pareillement au systeme GSM, le comportement du 
TEB a ete evalue pour des tranches temporelles 
composees de N=14 8 symboles chacune . Toutefois, les 
modulations retenues sont les modulations MDP4 et MDP8, 

15 Le canal discret considere est compose de un (L=l) ou 
de deux (L=2) ou de quatre (L=4) trajets de meme 
puissance moyenne. Deux configurations ont ete 
considerees : 

■ Une sequence de reference centrale de 26 symboles 
20 bipolaires prise parmi les 8 sequences utilisees 

dans le systeme GSM, Ces sequences sont des 
extensions cycliques des 8 sequences CAZAC de 
longueur 16 [1, 11]. 

■ Une sequence de reference centrale ({-1, +1, +1, 
25 +1, -1} pour deux trajets et {-1, +1, +1, +1,-1, 

+1, +1} pour quatre trajets de 4+L-l symboles 
bipolaires, obtenue par extension cyclique de 
1' unique (a un signe et un decalage pres) 
sequence CAZAC bipolaire de longueur 4 a savoir 
30 {-1, +1, +1, +1}. 
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Pour la figure 4, les modulations retenues sent les 
modulations MDP4 et MDP8 utilisant la premiere 
configuration avec deux trajets (L=2) o Les courbes se 
rapportent aux cas suivants : 

- 51 : MDP8, recepteur classique, 

- 52 : MDP8^ recepteur selon 1* invention, 

- 53 : MDP8, recepteur avec canal discret 

parf aitement connu, 

- 54 : MDP4, recepteur classique, 

- 55 : MDP4, recepteur selon 1' invention, 

- 56 : MDP4, recepteur avec canal discret 

parf aitement connu. 

Sur la figure 5, les modulations retenues sont les 
modulations MDP4 et MDP8 utilisant la deuxieme 
configuration et un canal discret avec deux trajets 
(L=2) • Les courbes se rapportent aux cas suivants : 

- 61 : MDP8, recepteur classique, 

- 62 : MDP8, recepteur de 1' invention, 

- 63 : MDP8, recepteur avec canal discret 

parf aitement connu, 

- 64 : MDP4, recepteur classique, 

- 65 : MDP4, recepteur de 1' invention, 

- 66 : MDP4, recepteur avec canal discret 

parf aitement connu • 

Pour la figure 6, la modulation retenue est la 
modulation MDP4 avec la deuxieme configuration et un 
canal discret avec quatre trajets (L=4). Les courbes se 
rapportent aux cas suivants : 
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- 71 : recepteur classique, 

- 72 : recepteur de 1' invention, 

- 73 : recepteur avec canal discret parf aitement 

connUo 

5 

Les courbes des figures 4 et 5 montrent que, a la 
difference d'un recepteur selon 1' invention, les 
performances du recepteur classique du systeme GSM sont 
tres sensibles a la longueur de la sequence de 

10 reference utilisee et se degradent substantiellement 
quand on passe de la premiere configuration a la 
deuxieme. Les courbes des figures 5 et 6 montrent, en 
outre, que les degradations en performances du 
recepteur classique augmentent signif icativement avec 

15 le nombre de trajets alors que, a 1' inverse, ces 
degradations restent presque inchangees et faibles pour 
un recepteur selon 1' invention, Le recepteur selon 
1' invention presente done une excellente robustesse vis 
a vis de toute reduction de la longueur de la sequence 

20 de reference et de toute augmentation du nombre de 
trajets du canal, 

Le tableau 2 caracterise de maniere chiffree, dans 
le cas de la modulation MDP4, 1' amelioration en 
25 performances apportee par le recepteur selon 
1' invention vis-a-vis de I'algorithme ameliore du 
systeme GSM. Le rapport BdTs est egal a 1/2500, le taux 
d'erreur binaire est egal a 10"^ et N=148. 
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Gain en dB 


Nombre de 
trajets ( Z, ) 


Sequence de 
r6f6rence 
longue 


Sequence de 
reference 
courte 


R6cepteur selon 1' invention avec 

s6quence courte, recepteur 
classique avec sequence longue 


1 


0,03 


1.03 


0,03 


2 


0,48 


1,83 


0,48 


4 


1.42 


3,04 


1.43 



Tableau 2 



5 Ce tableau montre que, pour une longueur donnee 

(premiere ou deuxieme configuration) de la sequence de 
reference, un recepteur selon 1' invention a toujours de 
meilleures performances qu'un recepteur classique 
ameliore. Ces resultats montrent aussi que 

10 1' amelioration apportee par 1' invention est d'autant 
plus grande que la longueur de la sequence de reference 
est reduite ou que le nombre de. trajets du canal est 
augmente. Ces resultats montrent enfin que, meme avec 
une sequence de reference courte (deuxieme 

15 configuration) , un recepteur selon 1' invention a 
toujours de meilleures performances qu'un recepteur 
classique avec une sequence de reference longue 
(premiere configuration) . 

20 La figure 7 represente le taux d'erreur binaire en 

fonction du rapport £5 / A^o pour un etalement Doppler 
normalise BdTs=1/100 contraignant et des tranches 
temporelles composees de N=100 symboles chacune. La 
modulation retenue est toujours la modulation MDP4. Le 

25 canal discret considere est compose de un (L=l) ou de 
deux (L=2) ou de quatre (L=4) trajets de meme puissance 
moyenne , Une seule configuration avec trois sequences 
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de reference reparties est consideree. Chaque sequence 
est composee de 4+L-l symboles bipolaires et est 
obtenue par extension cyclique (a un signe et un 
decalage pres) de 1' unique sequence CAZAC bipolaire de 
longueur 4 a savoir {-1, +1, +1, +1}. 

La valeur de I'etalement Doppler consideree est 
suf f isamment contraignante pour ne pas permettre a un 
recepteur classique de suivre les variations de chacun 
des trajets du canal discret qui peuvent subir des 
rotations de phase de ±7t entre le debut et la fin de 
chaque tranche temporelle . 

Les courbes de la figure 7 se rapportent aux cas 
suivants : 

- 81 : L=l, recepteur avec canal discret 

parf aitement connu, 

- 82 : L=l, recepteur de 1' invention, 

- 83 : L=2, recepteur avec canal discret 

parf aitement connu, 

- 84 : L=2, recepteur de 1' invention 

- 85 : L=4, recepteur avec canal discret 

parf aitement connu, 

- 86 : L=4, recepteur de 1' invention, 

Les courbes de la figure 7 raontrent que le 
recepteur selon 1' invention conduit a une degradation 
tres faible des performances vis-a-vis d' un recepteur 
ideal connaissant parf aitement le canal discret. Cette 
degradation est de I'ordre de 0,4 dB pour un seul 
trajet (L=l) , 0,5 dB pour deux trajets (L=2) et 1 dB 
pour quatre trajets (L=4) . 
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Les performances du procede de 1' invention ayant 
ete exposees, la phase de reception et notamment le 
traitement du signal regu, vont etre decrits plus en 
details . 

Un recepteur selon 1' invention determine^ 
conditionnellement aux echantillons regus du bloc a 
traiter, la realisation du canal discret la plus 
vraisemblable durant une tranche temporelle donnee. 
Cette determination se fait selon un critere dit de 
probabilite a posteriori maximale (MAP en abrege pour 
"Maximum A Posteriori probability" en anglais) . Cette 
estimation du canal necessite la connaissance des 
probabilites a priori des symboles des sequences de 
donnees transmiseso 

Par souci de simplification, les developpements 
mathematiques qui suivent et les figures 
correspondantes se rapportent toujours au cas ou chaque 
bloc a traiter pour 1' estimation de canal est compose 
uniquement des echantillons regus se rapportant a une 
seule tranche temporelle. Le cas ou les symboles (de 
reference ou non) d'autres tranches temporelles (du 
meme utilisateur ou d'autres utilisateurs/ du canal 
commun de signalisation, ...) sont utilises se traite 
de la meme fa?on. 

On note I'operateur de transposition. Comme 

indique en figure 2, toute tranche temporelle d'un 
utilisateur donne est composee d'une suite de N 
symboles aQ,ai,...,ajs4-i que 1 ' on peut ecrire sous forme de 
vecteur colonne 
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les symboles etant eventuellement modules en phase 
(MDP2, MDP4, MDP8, ,..). L'energie transmise 
associee k chaque symbole est ind^pendante des 

valeurs prises par ce symbole mais peut varier d'un 
symbole de la tranche temporelle a 1' autre. Cette 
variation est generalement la consequence d'une 
allocation par le systeme AMRT d' energies differentes 
aux sequences de donnees et aux sequences de reference. 
Comme indique egalement dans la meme figure, les 
echantillons en sortie du filtre blanchisseur 
correspondant a cette tranche temporelle forment un 
vecteur colonne 

de N-^-L-l composantes, ou L est le nombre de trajets 
du canal discret. 

La figure 8a rassemble les moyens compris entre 
1' entree du modulateur de I'emetteur et la sortie du 
filtre blanchisseur du recepteur. Cette figure fournit 
done une representation, du canal radioelectrique reel 
et du traitement effectue au niveau du recepteur 
(demodulation, filtrage adapte et filtrage 

blanchissant) . Sur cette figure, on trouve le 
modulateur 91, un bloc 92 symbolisant le canal multi- 
trajets, le demodulateur 93 avec . son filtre adapte 94 
et son filtre blanchisseur 95. 

La figure 8b montre une representation de ces 
moyens sous forme de registre a decalage avec (L-1) 
cellules Yi, Y2, Yl-i- Ce registre revolt le symbole 

courant et les cellules contiennent les symboles 

retardes d'une periode symbole, a savoir an-i, an-2/ ..-/ 
^n-(L'i)* Le circuit represents comprend encore (L) 
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multiplieurs Zq, 2i, Zl-i recevant respectivement 

(L) coefficients cj, , ... cl;'\ un additionneur S 

relie aux multiplieurs, un dernier additionneur 
recevant un signal Wn representant un bruit. Cet ultime 
additionneur delivre un signal rn- 

Chaque echantillon regu peut se mettre sous la 

forme 

L-\ 

A: = 0,l,...,A^+L-2, 

/=o 

ou cl est le / eme coefficient du filtre de canal 
discret au A: erne instant d' echantillonnage et est un 
bruit additif blanc gaussien. Pour simplifier la 
presentation et avoir des equations concises, les 
coefficients du canal discret, verifiant k<l ou 

k>N'^l, interviennent dans 1' expression de R^^ alors 
qu' ils sont nuls. Pour la meme raison, les symboles 
emis, verifiant k<Q ou k>N-\, apparaissent dans 
cette expression de R^^ meme s'ils sont nuls. 

Les coefficients du canal discret sont presque 
decorreles entre trajets discrets differents 
(independance vis-a-vis de la variable / ) . Abstraction 
faite des effets de la modulation et de la 
demodulation, cette propriete est due a 1' independance 
des trajets multiples du canal radioelectrique reel. 
Toutefois, les coefficients d'un meme trajet discret / 
sont correles entre eux (dependance vis-a-vis de la 
variable k) . Leur fonction d' autocorrelation discrete, 
definie par 
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depend uniquement du spectre de puissance Doppler du 
tra j et correspondant . 

Les spectres de puissance Doppler generalement 
rencontres en communications radiomobiles sont soit de 
type classique, (en environnements exterieurs), soit de 
type plat, (en environnements interieurs) . La fonction 
d' autocorrelation du /erne trajet, de puissance moyenne 
est donnee par 

dans le cas d'un spectre de puissance Doppler classique 
et par 

iPl = 4 siniTTkB^Ts ) / rrkB^Ts , 
dans le cas d'un spectre de puissance Doppler plat, ou 
Jo(') est la fonction de Bessel de premiere espece 
d'ordre 0 [3, 12]. La puissance moyenne varie d'un 
tra jet a 1' autre et caracterise le profil de 
I'intensite des trajets ("multipath intensity profile" 
en anglais [3]) du canal discret. 

Le recepteur de 1' invention doit determiner selon 
le critere du maximum a posteriori (MAP) , la 
realisation la plus vraisemblable des vecteurs 

<^'=(^/.^Li,^L..-.^iv.My / /= 0,1,.. .,1-1, 

caracterisant les L trajets discrets retenus par le 
recepteur pour 1' egalisation. 'A cette fin, il a besoin 
d'une representation adequate du canal discret lui 
permettant de s'affranchir de la correlation temporelle 
des coefficients d'un meme tra jet. II a egalement 
besoin d'une representation qui lui permette de 
s'affranchir de la variation eventuelle de I'energie 
transmise d'un symbole de la tranche temporelle a 
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1' autre. Cette representation est obtenue tout d'abord 
par 1' introduction du vecteur de symboles transmis 
normalises 

A = ( Aq , , , . . , Af^f_^ ) 

5 avec des composantes A^^ =^ I ^[E^ de module constant egal 
a un. Elle est egalement obtenue en reecrivant les 
echantillons regus durant une tranche temporelle sous 
la forme : 

10 ou C[ est la /:eme composante du vecteur normalise 

^ ~ (^/ V^o V^i >^/+2 V ^2 »--'>^A^+/-i yj^N'X y 
associe et caracterisant directement le /erne trajet du 
canal discret. Elle est enfin acquise grace a une 
decomposition, dite de FCARHUNEN-LOEVE [3] , de chaque 

15 vecteur normalise C' dans une base orthonormale (b*)^"' 

tenant compte a la fois de la variation de I'energie 
des symboles emis et des caracteristiques du spectre de 
puissance Doppler du /eme trajet [3] . Cette 
decomposition permet en effet de decrire le vecteur 
20 normalise C' (avec composantes temporelles correlees) 
du /erne trajet discret par un vecteur colonne. 
equivalent, G' = (Gg ,Gi',...,G!^.,)' , dont les N composantes 
sont bien decorrelees* Le vecteur peut etre 

reconstruit a partir de sa representation adequate G' 
25 par la transformation 
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Inversement, les composantes de cette 

representation adequate G' peuvent etre retrouvees a 
partir du meme vecteur C' par la transformation inverse 

Gi=(B*)-'C\ /: = 0,1,...,^-1, 
ou (•)* est I'operateur de conjugaison. 

Conformement aux proprietes de la decomposition de 
ECARHUNEN-LOEVE [3], les vecteurs de la base 
orthonormale associee au / eme trajet ne sont autres que 
les vecteurs propres de la matrice de covariance du 
vecteur normalise associe a ce trajet. De plus, les 
valeurs propres {^^^}^^^ associes ne sont autres que les 
variances des composantes Gl,G[,.,.,G^^_^ de la 
representation adequate G' du trajet correspondant . Les 
vecteurs {b*}^J associes a chaque trajet, ainsi que les 

valeurs propres correspondantes, peuvent done etre 
obtenus a partir de 1' equations matricielle 

H'B^=A*B*, A: = 0,1,...,A^-1. 
La (/w,«)eme entree de la matrice de covariance H' 

du vecteur normalise associe au /erne trajet est donnee 
directement par 

ipl^„ denote la fonction d' autocorrelation discrete de 
ce trajet definie plus haut . 

Le recepteur selon 1' invention doit normalement 
disposer de L bases orthonormales associees a chacun 
des L trajets du canal discret afin de reconstruire ou 
decomposer toute estimation intermediaire ou finale du 
canal discret. Dans la pratique, les formes des 
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spectres de puissance Doppler des L trajets du canal 
sont identiques. Les matrices de covariance Ih'}^"' 
associees sont done identiques, a un facteur 
multiplicatif pres, et les bases orthonormales 
correspondantes le . sont aussi. Une seule base 
orthonorinale est alors generalement suffisante au 
niveau du recepteur pour effectuer convenablement 
1' estimation du canal. 

Pour generer la base orthonormale du /erne trajet, 
le recepteur doit disposer d' une estimation aussi 
precise que possible de la matrice de covariance H' de 
ce trajet, A cette fin, il doit disposer d' une 
connaissance aussi precise que possible des 
caracteristiques du spectre de puissance Doppler du 
canal radioelectrique . Cette connaissance est mise a 
jour regulierement et comporte obligatoirement une 
estimation de I'etalement Doppler Bd ou d' une borne 
superieure de celui-ci. Elle peut comporter 
eventuellement une estimation de la puissance moyenne 
de chacun des trajets du canal discret et de la forme 
du spectre de puissance Doppler. Conditionnellement a 
cette connaissance, le recepteur peut, soit calculer en 
temps reel les bases orthonormales adequates, soit 
rechercher dans une banque de bases orthonormales pre- 
calculees celles qui concordent le mieux avec les 
caracteristiques du canal reel. En general, le 
recepteur n'a ni connaissance precise de la forme du 
spectre de puissance Doppler, ni connaissance parfaite 
et instantanee de I'etalement Doppler. Dans ce cas 
typique, il peut eventuellement faire I'hypothese d'un 
spectre de puissance Doppler plat avec, comme etalement 
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Doppler, une borne superieure de I'etalement Doppler 
reel. Par ailleurs, le recepteur peut ne pas disposer 
d'une estimation precise de la puissance moyenne de 
chacun des trajets du canal discret pris 
individuellement . Dans ce cas aussi, il peut aussi 
faire I'hypothese d'un canal discret avec L trajets de 
puissances moyennes egales, 

L' estimation des vecteurs c' caracterisant les L 
trajets du canal discret se fait iterativement selon le 
critere du MAP a travers les vecteurs normalises C' ou 
leur representation adequate G' . Cette estimation peut 
prendre en compte de maniere optimale aussi bien les 
caracteristiques des symboles de donnees et de 
reference emis (statistiques et energies) que les 
caracteristiques du canal multi-tra jets . 

L' estimation au MAP {g'}^^/J d'une realisation de la 

representation adequate {C'j^^J'J du canal discret est la 
valeur 

qui maximise la densite de probabilite conditionnelle a 
posteriori /7({g'}^^'J|R) . Conformement a 1' invention, et 
en utilisant 1' extension dite SAGE [7] de I'algorithme 
EM [8, 9], il est possible d'obtenir iterativement une 
solution aussi proche qu'il le faut de la ■ solution 
exacte . 

En general, la densite de probabilite 
conditionnelle a posteriori /?({g')^^JJ|R) a maximiser 
possede plusieurs maxima globaux conduisant souvent a 
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une ambiguite de 1' estimation de canal selon le critere 
MAP. Cette ambiguite peut etre levee grace a 
1' utilisation des symboles de reference connus du 
recepteur. Cependant, ceci se revele souvent 
insuffisant, car cette densite de probabilite 
conditionnelle possede egalement des maxima locaux 
pouvant etre atteints par I'algorithme SAGE a la place 
de 1' unique maximum global. Pour eviter ce probleme, il 
faut utiliser les symboles de reference pour determiner 

convenablement les conditions initiales {g'^**^}^^_"J , 

L'algorithme SAGE reestime successivement , par 
induction, chacun des vecteurs en gardant inchangees 
les estimations precedentes des autres vecteurs G" , 
avec m^l . Conformement aux proprietes de l'algorithme 
SAGE, ces reestimations successives garantissent une 
croissance monotone de la densite de probabilite 

conditionnelle a posteriori x{g')^'''J|R) . 

Etant donne le vecteur R d' echantillons re?us 
durant une tranche temporelle a traiter, le recepteur 

commence par calculer des conditions initiales (g'^*'^K'' 

des vecteurs {G'}^^JJ a partir des echantillons ne 
dependant que des sequences de reference. A cette fin, 
les sequences de reference doivent etre telles qu'une 
separation rapide et simple des contributions des L 
trajets a chacun de ces echantillons soit garantie. 
Pour pouvoir effectuer une estimation initiale, 
permettant le lancement de l'algorithme d' estimation de 
canal, il faut recourir a une ou plusieurs sequences de 
reference ayant chacune de bonnes proprietes 
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d' autocorrelation, A I'instar du systeme GSM, ces 
sequences peuvent etre construites a partir de 
sequences CAZAC par extension ou permutation 
circulaire. Ces sequences sont composees de symboles de 
module constant et possedent une fonction 
d' autocorrelation circulaire nulle partout a 
1' exception de I'origine. Dans le cadre de la presente 
invention, il est preferable d' avoir des sequences 
CAZAC avec des symboles appartenant a la constellation 
de la modulation utilisee pour vehiculer les symboles 
de donnees. A titre d'exemple, pour une modulation 
MDP2, la seule sequence CAZAC (a un signe et a une 
permutation circulaire pres) la plus courte, a savoir 
{ "1/ +1/ +1/ +1} / ne compte que quatre symboles- Pour une 
modulation MDP4, utilisant la constellation (l,/ -1,-/} , la 
seule sequence CAZAC (tou jours a un facteur 
multiplicatif i,-l ou et a un decalage circulaire 

pres) la plus courte, a savoir {l,/} , n'est composee que 
de 2 symboles. 

Une sequence de reference basee sur une sequence 
CAZAC de longueur M ne permet d'estimer que M trajets 
du canal discret au plus. Elle est obtenue par une 
extension cyclique de L-1 positions de cette sequence 
CAZAC et comporte au total M+L-1 symboles. La longueur 
M de la sequence CAZAC doit done etre superieure ou 
egale au nombre de trajets L du canal discret a 
estimer. Dans le cadre de la presente invention, il est 
parfois souhaitable d'utiliser des sequences CAZAC 
courtes pour nombre de trajets L a estimer petit 
(inferieur ou egal a 4 pour une modulation MDP2 et a 2 
pour des modulations MDP4, MDP8, -..), car elles 
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permettent de reduire signif icativement le nombre total 
de symboles de reference par tranche temporelle, tout 
en garantissant un echant illonnage plus regulier et 
plus frequent du canal discret, surtout dans les 
environnement a grands etalements Doppler. 

La figure 9 illustre 1' estimation d' un canal 
discret compose de 2 trajets {L=2) dans le cas d'une 
modulation MDP2 . La premiere ligne. A, donne la 
sequence de reference, la seconde, Tl, correspond au 
premier trajet, la troisieme, T2, au second trajet et 
la quatrieme, S, donne la composition des echantillons 
re?us. La sequence de reference elementaire utilisee 
{-1/ +lf +1, +lr -1}, composee de 5 symboles (M+L-l=5) 
n'est autre que 1' extension circulaire d'une seule 
position de la sequence CAZAC {-1, +1, +1, +1} composee 
de 4 symboles (M=4). Comme I'indique la figure 5, la 
sequence de reference elementaire occupe les positions 
entre n-4 et n de la tranche temporelle a traiter. Les 
echantillons regus, dependant uniquement de cette 
sequence, occupent quant a eux les positions entre n-3 
et n. lis sont donnes explicitement par : 

R„ = c'„A„ + c:,A„., +iv„= -c° +c:+w„ 

Pour une estimation peu biaisee des coefficients 
C° et Cl du canal discret normalise, pour k pris entre 
n-3 et n, il est necessaire de faire I'hypothese d'un 
etalement Doppler suffisamment faible pour garantir une 
variation negligeable de canal entre les positions n-3 
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et n , Dans ce cas, une estimation approximative de ces 
coefficients est obtenue en correlant les echantillons 
' ^ pris entre w-3 et n ^ avec la sequence CAZAC 
elle meme et un decalage circulaire d'une position vers 
la gauche de celle-ci. La correlation avec la sequence 
CAZAC {-1, +1, +1, +1} donne : 

comme premiere estimation des echantillons n-3,n-2, n-1 
et n du deuxieme trajet, ou est un terme qui tient 

compte du bruit et du biais. La correlation avec la 
sequence obtenue apres decalage circulaire, {+1, +1, 
+1, -1}/ donne quant a elle 

comme premiere, estimation des echantillons aux instants 
n-3,n-2,n-l et n du premier trajet, ou est un 

terme qui joue le meme role que . 

Dans le cas general, 1' estimation grossiere des L 
coefficients du canal, a un instant donne de la tranche 
temporelle, est obtenue par correlation des 
echantillons dependant uniquement de la sequence de 
reference avec L versions decalees de la sequence CAZAC 
correspondante , 

Chaque sequence de reference construite par 
extension circulaire d'une sequence CAZAC permet 
d'obtenir un echantillon bruite de chacun des trajets 
du canal discret- Le nombre d' echantillons par tranche 
temporelle necessaires a la reconstruction du canal 
depend de I'etalement Doppler normalise BqTs et du 
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nombre N de symboles composant cette tranche. Plus le 
produit BdTsN est grand, plus le canal varie d'une 
extremity a 1' autre de la tranche temporelle et plus le 
nombre de sequences de references elementaires doit 
5 etre grand. Une sequence de reference elementaire est 
un sous-ensemble de symboles groupes de la sequence de 
reference associee a la tranche temporelle obtenu a 
partir d'une extension circulaire d'une sequence CAZAC. 
Conformement aux proprietes de la decomposition de 

10 KARHUNEN-LOEVE du canal discret, plus une valeur propre 
A* est grande pour un trajet / donne^ plus la 
contribution Gj(B'* du vecteur propre associe a la 
representation adequate de ce trajet est importante en 
moyenne. La figure 10 montre la repartition des valeurs 

15 propres les plus importantes pour une tranche 
temporelle de N=100 symboles de meme energie normalisee 
a 1 et des etalements Doppler normalises BqTs egaux 
respectivement a 1/50 (courbe 113), 1/100 (courbe 112) 
et 1/500 (courbe 111) . Cette figure montre clairement 

20 qu'une representation convenable de chacun des trajets 
du canal discret peut s'appuyer, sans perte perceptible 
de performance, sur les quelques vecteurs propres de 
plus grandes valeurs propres. Pour des valeurs faibles 
du produit BdTsN, le vecteur propre de plus grande 

25 valeur propre suffit a lui seal a une representation 
precise de chacun des trajets du canal discret. Ce 
vecteur propre varie alors tres lentement au cours du 
temps et permet done, grace a une sequence de reference 
groupee au milieu de chaque tranche temporelle, de 

30 reconstruire un canal discret quasi constant sur toute 
la tranche temporelle. Ceci est le cas du systeme GSM 
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ou les etalements Doppler normalises sont supposes etre 
tres faibles at le canal est suppose etre constant sur 
une tranche temporelle donnee. 

Pour des etalements Doppler normalises un peu plus 
eleves, les deux vecteurs propres associes aux deux 
plus grandes valeurs propres sont necessaires a une 
reconstruction precise du canal discret. Ces deux 
vecteurs variant suffisarament pour qu'une sequence de 
reference composee de deux sequences elementaires , 
placees de part et d' autre du milieu de la tranche 
temporelle^ se revele necessaire, 

Dans le cas general, plus grand est le produit 
BqTsN, plus grand est le nombre de vecteurs propres a 
utiliser dans la representation adequate du canal 
discret. De plus, plus grand est ce produit, plus 
importantes sont les variations de ces vecteurs propres 
et done plus important est le nombre de sequences de 
reference elementaires a mettre en oeuvre. De surcroit, 
a longueur de sequence de reference constante, plus le 
nombre de sequences de reference elementaires est 
grand, plus la longueur de chacune doit etre diminuee. 
Pour un nombre de trajets L donne, la longueur d'une 
sequence elementaire, M+L-1, peut au mieux etre reduite 
a 2L-1 en choisissant une sequence CAZAC de longueur 
L ; d'ou I'utilite d' avoir des sequences CAZAC tres 
courtes, Selon 1' invention, la sequence CAZAC {-1, +1, 
+1, +1} de longueur M=4 est recommandee pour toutes les 
modulations de phase quand le nombre de trajets L est 
inferieur ou egal a 4 et que le produit BqTsN est 
faible. Selon 1' invention egalement, la sequence CAZAC 

{l,/} de longueur M=2 est conseillee pour un nombre de 
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trajets inferieur ou egal a 2 et des produits BqTsN 
assez importants- 

Soit K le nombre de sequences de reference 
elementaires par tranche temporelle a traiter. Solent 
Q'-^^+i,C/_^^2,-"3Q' ' / = 0,1,...,Z -1 , les echantlllons du canal 
discret procures par la /:eme sequence de reference 
elementaire et S^^ = {k - M -k-l.k - M ■¥2,...,k] I'ensemble des 
indices de ces echantlllons. Les conditions initiales, 
|q/(0)|^^J^ de la representation adequate du canal selon 

1' invention peuvent etre obtenues a partir de ces 
echantlllons en utilisant les equations : 

Le facteur de ponderation w^j' depend de la norme 
NI de la restriction des vecteurs propres aux seuls 
echantlllons en correspondance directe avec les 
echantlllons de canal procures par la sequence de 
reference. II depend egalement de la qualite de la 
contribution du vecteur propre mesuree par le 

rapport A!^ / , Plus la variance A'* du coefficient 
dans la representation du canal discret est grande par 
rapport a la variance Nq des echantlllons de bruits 
plus la ponderation correspondante est importante et 
s'approche de 1 . 

En se basant toujours sur le vecteur R 
d' echantlllons re^us durant une tranche temporelle, le 
recepteur selon 1' invention conduit ensuite une 
estimation iterative du canal discret en partant de ces 
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conditions initiales - Conf ormement a 

1' extension SAGE de 1' algorithme EM, I'estimateur de 
canal optimal selon le critere du MAP calcule 
iterativement, trajet apr^s trajet, la reestimation 

5 {g'^"*')}';; a partir de 1' estimation {g'^"'},';' . A chaque 
iteration, I'estimateur de canal calcule de proche en 
proche les vecteur G'^''*'^ / = 0,1,..., I - 1 , a partir des 
vecteurs {g et {g-^'^^JJ, calcules auparavant. 

Selon 1' invention, la ^eme composante de la (d + l)^me 

10 reestimation g'^'''"'^ du vecteur G' associe au /erne trajet 
et obtenue en se servant de 1' expression : 

^ /!=/ \ A, L /n=0 m^t-t-\ J 'J 

ou 

' 1 

et C^^'^^ est la weme composante de la reconstruction, 
C"'^'^\ selon la decomposition de KARHUNEN-LOEVE du meme 
trajet du canal discret a la d erne iteration- Cette 
reconstruction est donnee explicitement par 

Comme precedemment , les ponderations tiennent 

compte de la qualite de la contribution du vecteur B'* 
a la representation adequate du canal discret- Plus la 
contribution d' un vecteur propre est importante par 
25 rapport a la variance du bruit, plus sa ponderation est 
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proche de 1 et plus il joue un role important dans la 
reconstruction de 1' estimation du canal. 

Af in d' estimer la representation adequate, G'^''**^ , 
du /erne trajet en accord avec 1' expression precedente, 
5 le recepteur commence par calculer une estimation 

de C^-4„., en retranchant les contributions des autres 
trajets m^l \ Cette soustraction est faite dans 
I'hypothese ou le vecteur A„ = K.^,.,^,^„_^,.3j,...,/|J est 
10 emis, d'ou I'utilite de ponderer plus tard cette 
expression par la probabilite conditionnelle a 

posteriori ^(a„|r,{G"^^**>}^"J^,{g-^^>)^^^^ d' avoir reellement 

emis ce vecteur, Une fois cette estimation de A , 
realisee, le recepteur la multiplie par le conjugue A^,f 
15 de A„^i pour s'affranchir de la modulation. Cette 
operation est effectuee en supposant toujours que le 
vecteur A„ a ete emis, d'ou la necessite de ponderer 
1' estimation resultante : 

U-t c:^""' A„.„- z cr"' a) Ai, 

20 de par la probabilite conditionnelle a posteriori 

d' avoir emis ce vecteur, L' estimation finale 

K-I.c:^''''A„_„-f^cr^A„.VM^^ 

de est bien sur obtenue en prenant en compte la 

contribution de toutes les valeurs pouvant etre prises 
25 par A„ . 
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Une fois cette estimation , du «eme echantillon 
du /erne trajet obtenue pour n entre / est N-hl-l, il 
ne reste plus au recepteur qu'a la projeter sur chacun 
des vecteurs B* de la base orthonormale associee a la 
representation adequate de ce trajet. La projection sur 
le kem^ vecteur est ensuite ponderee par afin de 
tenir compte de la qualite de sa contribution a 
1' estimation du / eme trajet du canal a la (^/ + l)eme 
iteration. 

En pratique, seuls quelques vecteurs de base 
interviennent de maniere importante dans la 
representation de chacun des trajets. Dans ce cas, le 
nombre de vecteurs utilises dans la transformation de 
KARHUNEN-LOEVE et la transformation inverse peut etre 
reduit sensiblement sans perte perceptible de 
performances tout en garantissant une reduction 
importante de la complexite. 

L' estimation iterative du canal discret peut etre 
accomplie un nombre fini D de fois de telle maniere 

que 1' estimation globale {C^^^}^^^ (obtenue a partir de 
Jq^oJ^^J par la transformation de KARHUNEN-LOEVE) 
garantisse une degradation imperceptible des 
performances du recepteur par rapport a la solution 
optimale qui serait obtenue apres un nombre infini 
d' iterations . Cette estimation iterative suppose la 
connaissance des probabilites conditionnelles a 

posteriori ^(a„|r,{g"^^*»>};^^^^^^ pour toutes les 

valeurs possibles des vecteurs A„ ^ (A„,^,,,^,A„_^,,,^,..,,A„) , 
n = 0,l,.,.,A'' + Z-2 . Ces probabilites conditionnelles 
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peuvent etre calculees exactement par 1' algorithme de 
BAHL, connu aussi sous le nom d' algorithme BCJR [6]. 
Elles peuvent aussi etre calculees de maniere 
approximative et simple grace a des variantes 
linearises de I'algorithme de BAHL [10], 

A titre d' illustration, on presente ici les 
grandes lignes du calcul par I'algorithme de BAHL des 

probabilites conditionnelles ^(aJr,{G'"^^"'^}^"^q,{G"^'^^}^^^ . A 

cette fin, il faut remarquer, conformement a la figure 
8, que I'echantillon i?„ obtenu a 1' instant n depend S 
la fois du symbole A„ emis a cet instant et des L-1 
symboles ^„-i.^„-2.---»^„-(i-i) le precedant. Ainsi, 
1' ensemble modulateur-canal multi-tra jets-demodulateur 
peut etre modelise par une chaine de MARKOV a temps 
discreto Cette chaine de MARKOV est decrite par un 
treillis dont les etats representent a chaque instant n 
le contenu du registre a decalage de la figure 8 et 
dont les etiquettes des branches representent la 
version non bruitee des echantillons possibles a la 
sortie du filtre blanchisseur . A titre d'exemple, la 
figure 11 montre une representation en treillis d'un 
canal discret compose de L~2 trajets pour une 
modulation MDP2 a deux etats {-1, +1} et une longueur 
de tranche temporelle N=4 . 

Un tel treillis represente la sortie attendue du 
canal en fonction du temps et selon les symboles emis, 
Dans le cas illustre ou le canal est suppose presenter 
deux trajets, la sortie du canal est conditionnee par 
le symbole emis a 1 • instant considere et a 1' instant 
anterieur. Par exemple, en partant de 1' instant n=l^ on 
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peut emettre soit -1, soit +1, sachant qu'a 1 ' instant 
precedent 0, on n'a rien emis (noeud initial a 1 'extreme 
gauche). Si I'on a emis -1, on s'attend en sortie du 

canal a recevoir (-1) Co+(0)Co, soit -C^. quantite qui 
5 figure sur la premiere branche montante vers le noeud 
-1. Si I'on a emis +1, on s'attend a recevoir 
( + 1) Cq+0 Cq , soit Cq, quantite qui figure sur la 
premiere branche descendante vers le noeud +1. 

Pour n==l, si I'on emet -1 et si I'on avait emis -1 
10 a 1' instant anterieur n=0, on s'attend a recevoir 
(-1) C?+(-l)cl , soit -Cf-c} . Si I'on emet +1, on 
s'attend a recevoir ( + 1) C?-c{=C?-c[ . De meme, si I'on 
avait emis +1 a 1' instant anterieur, on s'attend a 
recevoir -C?+c{ si I'on emet -1, et c{*+cl si I'on emet 
15 +1. Ces quatre valeurs figurent sur les quatre branches 
a n=l- La tranche s'acheve par une emission de 0, et la 
quantite attendue est -C^ ou C4 selon qu*on a emis -1 
ou +1 a 1' instant anterieur n=3. Le treillis represente 
done bien, a chaque instant^ les quantites attendues en 
20 sortie du canal, selon les symboles envoyes. 

Conformement a une telle representation du canal 
discrete chaque branche du treillis determine une 
valeur unique du vecteur A„ = (^„.(^.,j,^„.(^_2j,...,^J , car, 

d'une part, les L-1 premieres composantes sont donnees 
25 par I'etat de depart S„ = K.^^.,^,^^.^,.^^,...,^^.,) de la 

branche et, d' autre part, les L-1 dernieres composantes 
sont donnees par I'etat de destination 

-2)»^n-(L-3)'---'^n) • ^® probleme de calculer les 
probabilites conditionnelles a posteriori 
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p(aJr,{g"<''*'^}^'J^,{G'"^^>}^;;) se ramene done au calcul des 

probabilites des branches du treillis 

conditionnellement au vecteur d' echantillons regus R 
et a 1' estimation intermediaire du canal discret 
normalise donnee par les vecteurs C"*^^**^ pour les / 
premiers trajets aw = 0,1,..,,/-! et les vecteurs C"^*'^ pour 
les l^'i derniers trajets /w = /,/ + 1 . A cette fin, 
les branches du treillis doivent etre etiquetees par 
les coefficients de cette estimation intermediaire et 
non par ceux du canal discret reelc 

L'algorithme de BAHL permet de calculer 
ef f icacement les probabilites conditionnelles a 
posteriori des branches du treillis en factorisant 
1' expression de la probabilite du vecteur regu 
conditionnellement a chacune des sequences A pouvant 
etre emises, afin de ne pas dupliquer toutes les 
operations communes a des sequences emises ayant le 
meme passe ou le meme avenir a un instant donne. A 
cette fin, cet algorithme fait appel a des fonctions 
<ar„(S) et /?,(S) qui donnent respectivement un apergu du 
passe et de I'avenir pour chaque instant n==0, 1, 
N+L-1 et chaque valeur S possible de I'etat S„ . 

La fonction a„(S) est la probabilite conjointe 
d'etre a I'etat S„ = S a 1' instant n et d' avoir re?u les 

echantillons • La factorisation des calculs pour 

des sequences ayant un passe commun permet de calculer 
cette fonction grace a la recurrence : 

a„(S)=Z«n-.(S')/„(S'.S) 



SP 16348. C/RS 



o^i rn(S',S) est la fonction qui donne la probabilite, 
conditionnellement a un etat de depart du treillis 
S„-i=S', d'aller a I'etat S„=S et de recevoir en meme 
temps 1' echantillon Cette expression montre que 

1' information sur le passe peut etre acquise a la fois 
a travers les probabilites associees a 1' instant 
precedant immediatement 1' instant actuel et a travers 
les probabilites de transitions des branches du 
treillis. 

La fonction J3„(S) a une signification equivalente a 
celle de a„(S) . C'est la probability d' avoir regu les 

echantillons \ conditionnellement a I'etat de 

depart S„=S a 1' instant n. Cette fonction peut etre 
aussi calculee de proche en proche a partir de I'avenir 
par la recurrence 

y9„(S) = i;/5...(S)x„.,(S,S'). 

s* 

La encore, cette recurrence prouve que la 
connaissance de I'avenir pour un etat S„ = S peut etre 
obtenue a partir des caracteristiques a la fois des 
etats S„^|=S' suivant immediatement cet etat et des 
branches du treillis les reliant ensemble. Pour une 
branche du treillis caracterisee par I'etat de depart 
S„ =S' = (^„.(^,,j,>J„.(^_2j,...,>i„.j) et I'etat d'arrivee 

S„+i =S = , la fonction x*(S',S) est donnee 

par : 



r„.,(S',S) = PK)-^exp 



f 



1 



/n=0 ma/ 
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ou P(A„) est la probabilite a priori d'envoyer le 
symbole A„ . Ces probabilites a priori peuvent 
caracteriser un codage eventual des symboles modules 
reellement emis, lis peuvent egalement caracteriser les 
symboles de reference qui sont connus du recepteur et 
qui prennent une seule valeur possible. Pour ces 
symboles de reference, la probabilite a priori P(A„) 
prend la valeur 1 pour la valeur reellement emise et 0 
pour les autres valeurs possibles de la modulation. . 

Afin d'effectuer les recurrences sur les fonctions 
a„(S) et A(S) r il faut prendre comme valeurs 
initiales : 

^o(So = (0,0,.. .,0)) = 1 et >S;v..-,(S;v*.-, = (0,0,.. .,0)) - 1 

A partir du calcul par I'algorithme de BAHL des 
fonctions a„(S) et J3„(S) , il est facile d'obtenir les 
probabilites des branches du treillis, 

conditionnellement aux echantillons regus et a 
1' estimation intermediaire du canal discret, grace a la 
formule : 

pIa = aIrIc"^^*'^}'"' f C"^^^ ]'''']- ^" (S\S)^,,, (S) 
^ " ^ ^--'Z A)(So=(0,0,...,0)) 

ou A = (^„-(t-l)i^rt_(/.-2)»---'^n) f S' = (-^n-(Z.-l)»^/i-(f.-2)»--->^n-l) 

S = (-4„_(£_|),...,-^rt-l»-^n) • 

II faut remarquer toutefois que la complexity de 
I'algorithme de BAHL reste tou jours assez elevee et 
qu'en pratique il vaut mieux recourir a des algorithmes 
approximatif s qui permettent de diminuer notablement la 
complexite des calculs sans trop degrader les 
performances. Le recours a ces algorithmes 
approximatif s est d'autant plus conseille que le nombre 
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de trajets ou/et d'etats de la modulation est 
important . 

Apres D iterations, le recepteur dispose d' une 
estimation du canal discret {C*^^^}^^^ proche de la valeur 
optimale selon le critere du maximum a posteriori qui 
est normalement obtenue apres un nombre infini 
d' iterations. A partir de cette estimation du canal et 
des echantillons regus durant la tranche temporelle a 
traiter, le recepteur peut fournir soit des decisions 
fermes grace a I'algorithme de VITERBI, soit des 
decisions ponderees grace a I'algorithme SOVA ou a 
I'algorithme de BAHL. 

L'algorithme de VITERBI cherche la sequence emise la 
plus vraisemblable. II permet done de reduire la 
probabilite d'erreurs de decision sur ces sequences 
emises. L'algorithme SOVA est un algorithme de VITERBI 
avec sorties ponderees- II permet une caracterisation 
plus precise de la fiabilite des decisions effectuees 
par I'algorithme de VITERBI. L' algorithme de BAHL 
permet de determiner la densite de probabilite a 
posteriori de chaque symbole emis, conditionnellement 
aux echantillons re?us et a 1' estimation de canal 
obtenue. II permet done de determiner la valeur la plus 
vraisemblable de chaque symbole emis et done de reduire 
le taux d'erreur, Les sorties ponderees sont obtenues a 
travers les probabilites des symboles 

conditionnellement aux echantillons regus. 

Dans tous les cas, les sorties ponderees ne sont 
utiles que si un decodage correcteur ou un codage de 
source avec entrees ponderees est utilise afin 
d' augmenter les performances du recepteur. 
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A titre d' illustration, dans le cadre de 
I'algorithme de BAHL, le recepteur utilise 1' estimation 

du canal discret pour calculer les fonctions 
a„(S) et >^„(S) . Ces fonctions permettent alors de 
5 determiner pour chaque symbole A„ sa distribution de 
probability par la formule : 

p{a. = .^R.{G'<«),v:) = E i-fA. = a|r.{g«'),';:) . 

Toutes les operations precedentes d' estimation 
iterative du canal discret et de detection sont 
10 illustrees sur la figure 12. L'estimateur de canal 120 

utilise des moyens 121 definissant les bases 
/ = 0,1,...,Z-1 , et des moyens 122 calculant les 
coefficients de ponderation ®^ {^*}*l~o ' 

L'estimateur 120 met en oeuvre une premiere estimation 
15 grossiere du canal a partir des seules sequences de 
reference elementaires par correlation circulaire avec 
les sequences CAZAC correspondantes • Cette premiere 
estimation est symbolisee par le bloc 3q reference 
1201, qui delivre la representation adequate de 

20 1' estimation initiale . Ensuite, I'estimateur 

120 met en oeuvre une reestimation d'ordre d, 
representee par le bloc 3^ reference 1202, qui delivre 
la representation adequate de 1' estimation de canal 
j^^^J ^ suivie d'une reestimation d'ordre d + l, 
25 representee par le bloc 3^^, reference 1203, qui delivre 

la representation de 1' estimation de canal • 
Enfin, I'estimateur 120 met en oeuvre une derniere 
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reestimation d'ordre D, representee par le bloc 3^^ 
reference 1204, qui delivre la representation adequate 
de 1' estimation de canal {G'^^^j^^J'J . A partir de cette 
representation adequate du canal estime, I'estimateur 
fournit par les moyens 1205 1 ' estimation definitive 

du canal discret en mettant en oeuvre la 
transformation de PCARHUNEN-LOEVE. 

La f ourniture a un decodeur correcteur d' erreur ou 
a un codeur de source de decisions fermes ou souples 
sur les symboles emis, s'effectue par des moyens 123 
qui mettent eventuellement en oeuvre I'algorithme de 
VITERBI [3] , I'algorithme SOVA [2] , I'algorithme de 
BAHL [6] ou I'algorithme de BAHL linearise [10]. 

Selon 1' invention, la reestimation iterative du 
canal est effectuee D+\ fois et permet de calculer 

{^'^°^}/»o ' {^'^'^}/to ' • • • ' (^'^^^l/to' ' successivement . La 
figure 13 illustre les L etapes composant chaque 
iteration et qui sont necessaires a 1' estimation 
successive de chacun des L trajets. Grace aux moyens 
121 et 122 de la figure 12, (fournissant les bases 

{^'"^ir^o ponderations {y^^V^Io ^ ^ la (rf + l)eme 

iteration 3d+i de I'estimateur de canal met en oeuvre 
une premiere etape, basee sur I'algorithme de BAHL et 
representee par le bloc BAHLq reference 130, pour 
delivrer les probabilites conditionnelles a posteriori 

/'(a„|r,{G'"^''^}^;^) a partir de la d erne reestimation du 

canal. A partir de ces probabilites conditionnelles, la 
(c/+l)eme iteration met ensuite en oeuvre une deuxieme 
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etape, basee sur 1' algorithms SAGE et representee par 
le bloc SAGEq reference 131, pour permettre la (d -hi) erne 

reestimation G^^"""*^ du premier trajet. 

A partir de la (c? + l)dme reestimation des /-I 
5 premiers trajets, la (d+l)eme iteration met en oeuvre 
une troisieme etape, basee de nouveau sur I'algorithme 
de BAHL et representee par le bloc BAHL, reference 132, 
pour delivrer les probabilites conditionnelles a 



10 reestimation des L-l derniers trajets du canal et de 
la (<^-i-l)eme reestimation des / premiers trajets de ce 
canal. Par la suite, la (d-i-l)eme iteration met en 
oeuvre une quatrieme etape, basee sur I'algorithme SAGE 
et representee par le bloc SAGE; reference 133, pour 

15 permettre la (£/ + l)eme reestimation G'^"'"'^ du / eme 
trajet , 

Enfin, a partir de la (^/ + l)eme reestimation des 
I-l premiers trajets, la (J+l)eme iteration met en 
oeuvre une derniere etape, basee sur I'algorithme de 
20 BAHL et representee par le bloc BAHL^., reference 134, 
pour delivrer les probabilites conditionnelles a 



posteriori />(aJr,{G'"^^*'^}^;|^,G^'-*>-<^M a partir de la d erne 



reestimation du dernier trajet du canal et de la 
((i + l)eme reestimation des L-l premiers trajets de ce 
25 canal. Par la suite, la (^ + l)eme iteration met en 
oeuvre une derniere etape, basee sur I'algorithme SAGE 
et representee par le bloc SAGE^,, reference 135, pour 



posteriori 




partir de la d erne 
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permettre la (i/-hl)eme reestimation g^^"'^'^''**^ du ieme 
tra j et . 

La premiere iteration, symbolisee par le bloc 3^ 
reference 1201 de la figure 12, ne peut pas mettre en 
oeuvre les algorithmes de BAHL et SAGE pour calculer la 
premiere estimation du canal discret. Selon 
1' invention, on utilise alors les symboles de reference 
dont les probabilites conditionnelles a priori et a 
posteriori sont bien connues du recepteur pour 
effectuer cette premiere estimation. Cette premiere 
estimation s'effectue, en correlant les echantillons 
regus, dependant uniquement de ces symboles de 
reference, avec les sequences CAZAC correspondantes . 
Elle met en oeuvre les moyens 122 de la figure 12 pour 

calculer les ponderations • 
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REVENDICATIONS 



1- Procede de communications radiomobiles du type 
a acces multiple a repartition en temps, dans lequel on 
alloue a plusieurs utilisateurs des tranches 
temporelles pour 1' emission et la reception de signaux 
radioelectriques, et dans lequel : 

- a 1* emission, on constitue, dans chaque tranche 
temporelle allouee, des symboles de donnees et 
des symboles de reference, on emet un signal 
radioelectrique contenant ces symboles de 
donnees et ces symboles de reference, ce signal 
empruntant un canal radiomobile mult i-tra jets, 

- a la reception, on effectue un filtrage adapte 
du signal regu, on echantillonne le signal 
filtre au rythme des symboles, on effectue une 
estimation du canal radiomobile, on traite les 
echantillons en tenant compte de cette 
estimation du canal et I'on restitue les 
symboles de donnees propres a la tranche 
allouee, 

ce procede etant caracterise en ce que : 

- a 1' emission, pour chaque tranche temporelle 
allouee, on constitue plusieurs sequences de 
reference possedant des proprietes de correlation 
temporelle appropriees et I'on repartit ces 
sequences dans chaque tranche temporelle, 

- a la reception : 

i) pour estimer le canal radiomobile, on 
represente ce canal par un treillis 
comprenant des branches symbolisant les 
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transitions entre deux etats successifs du 
canal et caracterisant les sequences 
possibles a la sortie du canal, 
ii) on calcule les probabilites de ces branches 
en mettant en oeuvre un algorithme iteratif 
dans lequel on calcule, par un nombre fini 
(D+1) d* iterations et selon le critere de 
probabilite maximum a posteriori/ une 

estimation d'une representation {G^^^^j^Jd des 
L trajets de ce canal, 1' iteration de rang 
(d+1) (3d+l) permettant d'obtenir une 

reestimation ^^^"^^^^felo ^ partir de 
1* estimation {g'^*'^}^^ obtenue par 1' iteration 
precedente (3d) , et, a chaque iteration 
(3d+l), on calcule en L etapes les vecteurs 
Gf«i+i)^ ^=0, 1, L-1 a partir des 

vecteurs {g^^^^^^^ et {g-^^^^ calcules a 

I'etape precedente, 1' estimation (c'^^^I^Jo des 
L trajets du canal etant finalement obtenue 
a partir de 1' estimation de sa 

representation {g'^^^}^^q . 

2. Precede selon la revendication 1, dans lequel, 
a chaque etape de chaque iteration, on met en oeuvre 
successivement : 

Dun algorithme de BAHL ou une de ses variantes 
pour delivrer des probabilites 

conditionnelles a partir de la ji^me 
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reestimation des L-^ derniers trajets et de 
la (d+l)^*'"® reestimation des i premiers 
trajets, 

ii) un algorithme SAGE pour delivrer une (d+1)^^"^^ 
reestimation du l'^^^^ trajet. 



3. Procede selon la revendication 1, dans lequel, 

ayant obtenu une estimation du canal discret {c'^^^}^=o 
et, compte tenu des echantillons regus durant la 
10 tranche temporelle traitee, on prend une decision sur 
la sequence des symboles emis en utilisant un 
algorithme optimal de BAHL, ou une de ses variantes, 
pour obtenir la probabilite a posteriori de chaque 
symbole emis. 



15 



20 



4. Procede selon la revendication 1, dans lequel 
on repartit K sequences de reference dans chaque 
tranche allouee avec K compris entre 2 (inclus) et 12 
(inclus) . 

5. Procede selon la revendication 1, dans lequel 
chaque sequence comprend un nombre de symboles de 
reference inferieur a 8, 



25 6. Procede selon I'une quelconque des 

revendications 1 a 5, dans lequel les sequences de 
reference possedent une fonction d ' autocorrelation 
circulaire nulle partout sauf a I'origine. 

30 7. Procede selon la revendication 6, dans lequel 

on utilise une sequence de reference {-1, +1, +1, +1} 
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ou une sequence derivee de celle-ci par permutation 
circulaire, ou changement de signe pour un nombre de 
trajets (L) inferieur ou egal a 4. 

5 8. Precede selon la revendication 6, dans lequel 

on utilise une sequence de reference {1,1} ou une 
sequence derivee de celle-ci par permutation 
circulaire, conjuguaison ou changement de signe, pour 
un nombre de trajets (L) inferieur ou egal a 2. 

10 
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